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Abstract 



We propose a new and efficient method of reconstructing quark mass ma- 
trices from their eigenvalues and a complete set of mixing observables. By 
a combination of the principle of NNI bases which are known to cover the 
general case, and of the polar decomposition theorem that allows to convert 
arbitrary nonsingular matrices to triangular form, we achieve a parameter- 
ization where the remaining freedom is reduced to one complex parameter. 
While this parameter runs through the domain bounded by the circle with 
radius R = \/ (mf — m^) / (mf — m^) around the origin in the complex plane 
one sweeps the space of all mass matrices compatible with the given set of 
data. 
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I. INTRODUCTION 



The four three-generation mass sectors of quarks and leptons belong to the deepest enig- 
mas of the standard model of strong and electroweak interactions. While there is a great 
amount of experimental information of steadily increasing accuracy, theoretical models of 
mass matrices and mixing matrices are scarce. In this unsatisfactory situation it seems of 
utmost importance to parameterize the available data in such a way that possible textures 
in the mass matrices of a given charge sector become visible in an unambiguous manner. For 
instance, in the case of quarks, there are 10 data, viz. the masses of the three charge +2/3 
quarks, the masses of the three charge —1/3 quarks, and four observables in the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, to be compared to 12 physically significant parameters 
in the mass matrices and M^ d \ Therefore, reconstructing mass matrices from the 

data really amounts to finding an optimal parameterization that exhibits the remaining 
two-parameter freedom in a transparent way. 

An important step in this direction was taken by Branco et al. who realized that the 
class of so-called nearest-neighbour-interaction (NNI)-bases for chiral states are economical 
but still completely general |1|, and may hence be used in attempts to reconstruct the mass 
matrices from the observed mass eigenvalues and the empirical mixing matrices. These 
authors also gave an explicit procedure for constructing mass matrices in an NNI basis, 
for arbitrary mass sectors of quarks. Unfortunately, their analysis involves solving cubic 
equations. Although soluble in principle, these equations are too cumbersome to solve and 
do not allow for a practical and efficient reconstruction. 

In this paper we propose a new method of reconstruction that avoids these shortcomings. 
We conjecture that this method and the appropriate parameterization are optimal in the 
sense of concentrating the remaining freedom in a single complex parameter whose domain 
of variation can be restricted to the interior of a circle in the complex plane. In particular, 
we succeed in reconstructing the mass matrices proper, up to unobservable changes of basis. 
This goes beyond, say, the work of Harayama and Okamura [0] who express the CKM matrix 
in terms of six parameters, with a two-parameter freedom. The way from their result to the 
mass matrices seems involved and not well suited for a practical analysis. 

We make use of the polar decomposition theorem for nonsingular matrices || 

M = TW , 

where T is a lower-triangular matrix and W is a unitary matrix which, when applying this 
formula to three generations of chiral quarks, can be absorbed in the right-chiral fields. If, 
in addition, we work in the class of NNI bases in which the (21)-element of T is seen to 
vanish, we still cover the most general case but get rid of all redundant quantities. More 
precisely, if 

ft{q) = M {q) M { ^ = f {q) f (9)t , q = u,d, 

are the "squared" hermitean mass matrices and = diag (m\ , m%, rnf) with mi = m u or 
nid etc. their diagonal forms, then, in any NNI basis, 

#(«) = WD^U, H {d) = UW C KMD {d) V ] CKU U . 
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The matrix U which is known analytically j|, depends on two complex parameters, say a 
and b defined in eq. (^6]) below. These parameters are related through a quadratic equation 
whose coefficients are elements of the matrix 

i.e. of a matrix that is obtained solely from experimental data! Solving for one or the 
other of them, say b = b(a), reduces the set of admissible mass matrices to the expected 
two-parameter freedom in the variable a. Progress achieved in this way is twofold: On the 
one hand, parameterization in terms of, say, a is analytically simple and transparent. On 
the other hand, the domain of variation of a and the formulae for the elements of and 
are such that the space of admissible mass matrices can be studied graphically and 
numerically, as a sweeps through all allowed values, in a quantitatively reliable manner. 
Although we have not done this yet, one can even follow the propagation of the error bars 
of the experimental input, in not too involved a procedure. 

The paper is organised as follows: In sect. 2 we review the choice of NNI bases and 
recapitulate the relevance of the polar decomposition theorem for the problem at stake. 
Sect. 3 which is the main body of our work, describes the explicit construction of the matrix 
U as well as its parameterization in terms of a, b, and the squared masses of the up-sector. 
The NNI condition is encoded in the quadratic equation ( |3"3"D below. In sect. 4 we discuss 
symmetries helpful in solving that constraint, and consider some limiting cases in order to 
illustrate the method. We also propose an expansion of the solutions in terms of a parameter 
which, due to the hierarchy in the quark masses, is numerically very small. The final sect. 5 
gives two examples of how a conjectured texture of mass matrices can be checked against the 
data in a simple and transparent manner, by converting the mass matrices to the (general) 
form studied here. It ends with a few conclusions. 



II. GENERAL NEAREST-NEIGHBOUR BASES AND TRIANGULAR 

MATRICES 

In the case of three generations the mass matrices and in the so called nearest 
neighbour interaction (NNI) basis are characterized by the following generic structure, 

/0 * 0\ 

M {u \ M (d) = * ★ . (1) 
\0 * ★/ 

Here the "^"-entries appearing on the r.h.s. of ([I]) are arbitrary non-vanishing, complex 
numbers. The interpretation of this particular choice of the mass matrices that is usually 
given by its proponents is the following: the (11)- and (22)-elements are put equal to zero, 
while letting the (33)-element be different from zero, with the idea of describing an initial, 
no-interaction situation where two quarks are massless and only one is massive. Furthermore, 
only neighbouring generations are allowed to interact, by assuming nonvanishing (12)-, (21)-, 
(23)-, and (32)-elements, but vanishing (13)- and (31)-elements. 

However, it has been known for a long time that this setting, although very tempting 
and intuitive at first sight, is ill-defined unless it is supplemented by further assumptions. 
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Indeed, Branco, Lavoura, and Mota showed in |l] that any set of admissible mass matrices 
{MW.M^J of the standard model, i.e. any set of two completely arbitrary non-singular 
3 x 3-matrices, can be transformed to the form given in (Jl]) without the need for any further 
assumption. In other words, in the framework of the minimal standard model where only 
left-chiral fermion fields participate in charged current weak interactions, the mass matrices 
in ([!]) are still completely general, the specific form (|l|) reflecting no more than a specific 
choice of chiral basis. This fact is mainly due to the observation that weak interactions of 
quarks remain unchanged under the simultaneous transformations 

-> M {u) = U^M (u) V u and 
M (d) -> M (d) = U^M {d) V d , (2) 

of the mass matrices, where U, V u and Vd are arbitrary unitary 3 x 3-matrices: On the 
one hand, the unitaries V u and Vd act on the right-chiral fields (uR,CR,tn) and (cIr, SR,bn), 
respectively, and, hence, can be absorbed by a redefinition of these unobservable fields, 
without loss of generality. On the other hand, the common unitary matrix U in (0) which acts 
on the left-chiral quark fields, drops out when calculating the physical Cabibbo-Kobayashi- 
Maskawa (CKM) matrix 

^CKM = U^U™ . (3) 

In (§ and uf are the unitary matrices which diagonalize the "squared", hermitean 
mass matrices 

#(<?) = M^M^t , q = u,d, (4) 
of the up- and down-sector, respectively, viz. 

uWffMuW = dmg(m 2 u ,m 2 c ,m 2 ) = , 

XjWhMxjW = diagK, ml m\) = £>W> . (5) 

Thus, as proved in [|T|, the structure (|1|) of the mass matrices corresponds to no more than 
a special, physically admissible choice of the electroweak basis. The interpretation sketched 
above fully rests on this choice and will no longer be valid in other electroweak bases. Unless 
this special choice is singled out by additional arguments that could stem, e.g., from physics 
beyond the standard model, the above interpretation loses its physical significance because 
physics, of course, must not depend on the choice of basis. 

There is an alternative derivation of the same conclusion which, at the same time, helps 
us to fix notations for our subsequent calculations. As was shown in |4j], 0, when dealing 
with questions of reconstructing mass matrices from the experimental data, i.e. from four 
independent absolute values of CKM matrix elements and the six quark masses, a most 
economic but nevertheless physically completely general parameterization of the quark mass 
matrices is given by triangular mass matrices, 



/ * 

T (u) , T (d) = I ★ ★ 



(6) 



This parameterization results upon exploiting the polar decomposition theorem for non- 
singular, but otherwise arbitrary matrices: 

M (q) = T (q) w ( 9 ) ^ q = U: d, with W {g) unitary. (7) 

Again, the unitaries W^ q \ q — u,d, can be absorbed by a redefinition of the right-chiral 
quark fields, without loss of generality. 

At the level of the squared mass matrices the NNI structure ([!]) of the mass matrices is 
equivalent to a vanishing (12)-element (and - because is hermitean - also a vanishing 
(21)-element): 

H[f = = H { 2 f , q = u,d. (8) 

Because of = M^M^ = T^T^t, q — u,d, a straightforward calculation shows that 
within the triangular parameterization the NNI condition @ reads 

fg) = 0, q = u,d . (9) 

Mimicking for a moment the (invalid) interpretation mentioned in connection with the NNI 
structure ([!]), eq. (j|) would suggest that there is no direct interaction between the first 
and the second generation while, in the basis ([1]), these evidently do mix, and, in fact, mix 
strongly. Direct interactions seem to be present between the second and third generations 
as well as between the first and third generations only, due to the non-vanishing (32)- and 
(31)-elements of T^ q \ respectively, in contrast to ([]]) where seemingly there is no direct cou- 
pling between first and third generations. 

These statements underpin once more that such an interpretation is dependent on the elec- 
troweak basis chosen for the representation of the mass matrices and, hence, should better 
be avoided altogether. 

Although it is unrelated to a specific physical picture of quark masses and mixings, bases 
that yield the NNI form of the mass matrices turn out to be very useful in the process of 
reconstructing mass matrices from the observed quark mixings and masses. Therefore, in 
what follows we shall make extensive use of this class of bases. That is to say, we start from 
triangular mass matrices whose (21)-elements are zero, viz. 

fa \ fa' \ 

f {u) = o Jo, f (d) =(o $' o . (10) 

Here and in the sequel the hat on the symbols refers to the choice of an NNI basis. 

Possible phases can be absorbed into the right-chiral fields and, hence, without loss of 
generality, the matrix elements a, (3, 7 and a', (3', 7' along the diagonals may be chosen to 
be real. In fact, also the phases (pi and <p 3 could be dropped by making use of the unitary 
matrix U in ([J) thereby rendering real. In doing so, on the theoretical side we are left 
with 12 parameters, 5 from the up-sector and 7 from the down-sector. These have to be 
confronted with 10 experimental data, i.e. 6 quark masses + 4 real observables in the CKM 
matrix, leaving a freedom of two parameters. This is characteristic for an NNI basis. 
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III. AN EFFICIENT RECONSTRUCTION PROCEDURE 



The authors of |1J] give a detailed prescription for the construction of mass matrices in 
an NNI basis, for arbitrary mass matrices: The main step consists in solving the eigenvalue 
problem 

(Hto + KHW)ij u j2 = \u i2 (11) 

for the matrix + nH^ d \ where k denotes an arbitrary complex number. Note that 
k reflects the two-parameter freedom within the class of NNI matrices. Once the second 
column (-Uj 2 ) of the unitary matrix U, which transforms the arbitrary squared mass matrices 
(|) to NNI form (see (H)) is determined according to (|ITD, the first column (un) is calculated 
by means of 

u {1 oc e ijk u* j2 u* l2 (H {u) )i k , (12) 

see |1[. Finally, the third column (wj 3 ) follows from the unitarity of U. It is easily checked 
that (HD and (0) indeed imply u* (# (9 % u j2 = for q = u, d as required. 



However, the prescription just outlined is not very well suited when aiming at an explicit 
construction of all NNI mass matrices. This is simply due to the fact that in the interesting 
case of three generations the eigenvalue problem ( |TTD leads to a cubic equation. The solutions 
of this cubic equation are, of course, known in principle but the corresponding expressions 
are rather lengthy and involved and they complicate tremendously subsequent calculations. 
In this paper we propose a different procedure which will be seen to result in much simpler 
expressions. 

We start from the following observation || : Without loss of generality we may assume that 
the squared mass matrix of w-quarks is already diagonal, that is, in other words, that the 
mixing has been shifted entirely to the down-sector. Indeed, this is achieved by exploiting 
once more the freedom contained in (|2]) by choosing U = . With this choice we obtain 
(using (§, ©, and (§): 

#(«) U^H^U™ = and (13) 

ff(tt) ! y fl{d)jj{u)\ _ jj{u)jj{d)\ ^d)jj{d)jj(u)\ 

L L L L L L 

= Vckm.D^V ckm 

Thus, in what follows, by this redefinition, the experimental input will be coded in the form 

HM = and = VckmD^V^ . (14) 

Note that is completely fixed in terms of the experimental input, i.e. the three quark 
masses of the down-sector and four physically relevant parameters of the CKM matrix. Let 
us comment on this point in more detail: Because, in general, normed eigenvectors are only 
fixed up to arbitrary phase factors, instead of W = we can as well choose = P\U^/ 
in the above reasoning, with Pi a diagonal phase matrix, 

P x = diag^e^ 2 ,^ 3 ) . (15) 
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With this choice the second substitution in fll3f) becomes: 

H {d) ^ ft(d> = Pi y cKM D (d) y£ KM pl 

= PiV CKM P ] 2 P 2 V ] ckm P\ (16) 

In the last line, P2 again denotes a diagonal matrix containing only phase factors, 

P 2 = diag(e^,e^,l) , (17) 

and the diagonal character of has been used. Eq. ([16]) shows that we are allowed to 
choose any parameterization for Vckm that we would like to start with. The transition to 
any other parameterization of Vckm is easily accomplished by means of suitable choices of 
the phase matrices Pi and P 2 . 

Next, we have to tackle the task of finding the unitary matrices U which transform the 
squared mass matrices to NNI form. Postponing for a moment the analysis of the 
down-sector, the corresponding condition for U in the up-sector simply reads 

ft(u) = rf D V) V D (u) = JJJjHjj^ _ ( lg ) 

In other words, at this stage U is given by those unitary matrices that diagonalise the 
most general squared mass matrix 

H (u) = T (u) T (u)] where T (u) jg taken from ^ ( with 

02 = = v?3, without loss of generality). This is the first condition for the matrix U. 

An analytical expression for the matrix U is obtained by restricting to the case k% = the 
general solution of the problem of diagonalisation that we had obtained earlier in [|J, viz. 

(f(m u )/ni g{m u )/ni h(m u )/ni\ 
f(m c )/n 2 g(m c )/n 2 h(m c )/n 2 , (19) 
f (m t )/n 3 g(m t )/n 3 h(m t )/n 3 J 

where the functions /(wj), g(mi), h(mi) and the denominators are given by 

f( m ) = -ak 3 2 -mf) "I 

g{rrii) = -(3k 2 (a 2 - mf) > (i = u,c,t) 

Kmi) = (a 2 - m 2 )0 2 - mf) J 

n\ = (a 2 - m 2 u )0 2 - m 2 u )(m 2 - m 2 u )(m 2 c - m 2 J 
n 2 = (ml — a 2 )0 2 — m 2 ){m 2 — m 2 ){m 2 — m 2 u ) 
n l = (m\ — a 2 ){m 2 — (3 2 )(m 2 — m 2 c ){m 2 — m 2 u ) . 

The diagonal phase matrix P, 

P = diagt-l^ 2 ^ 3 ) , (20) 

represents the freedom of multiplying each eigenvector of H^ u > with an arbitrary phase 
factor[|. As we are aiming at all NNI mass matrices this freedom must be taken into account. 



The choice of the first phase, exp{i7r}, is made for the sake of convenience. 
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This will become clear also in a moment when we count the degrees of freedom explicitly. 
Furthermore, from the comparison of the characteristic polynomials of T^T^ and of 
the parameters k 2 and k 3 are fixed in terms of a, (3 and the squared quark masses of the 
up-sector by means of the relations: 

m 2 u m 2 c m\ = a 2 (3 2 f (21) 
m 2 u m 2 c + m 2 u m 2 + m 2 c m 2 = a 2 (3 2 + (3 2 ^ 2 + tfa 2 + a 2 k\ + (5 2 k\ (22) 

m l + m l + m t = « 2 + P 2 + 7 2 + k \ + k l ( 23 ) 

As a consequence, U is a function of four real parameters (and, of course, the masses of the 
up-quarks) , 

U = U(aJ,ij 2 ,ij 3 ) , (24) 
see appendix A for more details. 

Next we turn to the down-sector. In order to guarantee the NNI form in the down-sector, too, 
the unitary matrix U, eq. fll9|), has to fulfill one additional condition. With = WH^U 
and setting U = (uy) and = VckmD^V^ = (^), see G14T), this condition reads 

(H&) 12 = u\ x h %] u 32 = . (25) 

This yields the second condition for the matrix U. As this is an equation for complex num- 
bers, two out of the four parameters a,$,ip2 and ip 3 are fixed this way, leaving a freedom of 
two real parameters. This is characteristic for the NNI form. Please note that it is essential 
to take into account properly the freedom parametrized in P, eq. ( p0|) . Had we missed the 
phase matrix P, U would have been completely determined by ( plf) and only one special 
set of NNI mass matrices would have resulted, contrary to our purpose of reconstructing all 
NNI mass matrices. 

The parameterization of U, eq. (|T9|), in terms of a, /3, ip2 and ip3 is not suited for constructing 
the solutions of (p5|), simply because the resulting equation contains a complicated sum of 
cosines with arguments ip 2 , as well as their difference ip2 — ' l p3- The practical reconstruc- 
tion in this framework would not be simpler than within the original proposal of Branco, 
Lavoura and Mota. The situation changes decisively if we use a new parameterization of U 
in terms of two complex numbers defined as follows^ 

a := and b := . (26) 

In particular, the moduli and phases of a and b are given by 

(27) 



'0 2 


-m 2 u ){m 2 


— a 2 )(m 2 




(a 2 


-m 2 u )0 2 


-m? c ){m 2 t 


-m 2 c ) 


'0 2 


-m 2 u )(m 2 


— a 2 )(m 2 


-m 2 ) 



1^1 _ / \p - " L u)\:" L t - " A'» e ~ in u) ^8) 

(a 2 — ml)(m 2 — (3 2 ){m 2 — m 2 ) 



2 The case u\ 2 = must be excluded at this point. However, as we shall see below, this is no 
restriction: The case where a tends to infinity, a — > oo, is mapped, by a symmetry of the equations, 
to the point a = 0, cf. eq. (^) below. 



S 




for m 2 < a 2 < ml < (3 2 < 
for m 2 u < (3 2 < ml < a 2 < m 



ml 



.2 ' 



(29) 
(30) 



The two complex variables a and b replace the four real variables a,/3,ifj 2 and ip 3 . In fact, 
a straightforward calculation shows that U, eq. (ITS), when expressed in terms of the new 
variables a and b reads as follows: 



U 



/ 


(m 2 -m=)|o||6| 


i 


a 2 (m 2 — m 2 ) + 6 2 (m 2 








iV 2 


iVi 










(m* -m*)|&|a 


a 


a(\b\ 2 (ml — ml) — (ml - 








iV 2 |a| 


iVi 














b(\a\ 2 (ml — ml) + (ml — 






V 


iv 2 |&l 


iVi 


AW 




) 



(31) 



with N 2 = 1 + |a| 2 + |6| 2 (32) 
and N 2 = \ab\ 2 (m 2 t - m\) 2 + \a\ 2 (m 2 c - m 2 ) 2 + \b\ 2 (m 2 t - m 2 ) 2 

As long as only the up-sector is under consideration a and b remain arbitrary and are not 
restricted at all. It is the second NNI condition (|25|) that imposes a constraint on them: in 
the new parameterization this condition takes the simple form of a quadratic equation, viz. 

(m\ — m 2 ) ab (hn + ah\2 + 6^13) 
- (m\ - m 2 ) b (h 2 \ + ah 2 2 + bh 2 z) 

+ (m 2 c - m 2 ) a (h 31 + ah 32 + bh 33 ) = (33) 

Depending on whether (|3"3] ) is solved for a = a(b) or for b = b(a) the complex parameter b 
or the complex parameter a remains free and, thus, we recover the two-parameter freedom 
of the NNI reconstruction. 

By means of the above formulae elementary calculations yield all parameters of the NNI 
mass matrices in terms of a and b. For instance, for (3 we obtain: 



m 2 u + 


a 


2 m 2 c + \b\ 2 ml 


1 + 


a 


2 + \b\ 2 



(34) 



After insertion of a = a(b) or b = b(a) according to ( pB]) this equation specifies all admissible 
values for the parameter (5 in the NNI form of the mass matrices by varying the unconstrained 
parameter b or a, respectively. 

The results for all other parameters are quoted in appendix B. 



IV. PARAMETER DEPENDENCIES, SYMMETRIES AND EXPANSIONS 

The results derived in the previous section whose details are spelled out in appendix B, 
are completely general and analytical in the sense that no approximations whatsoever have 
been made. Furthermore, we conjecture that the parameterization and reduction to the 
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complex parameter a (or, alternatively, the parameter b) described in the previous section 
is the best one can do in reconstructing mass matrices from their eigenvalues and the CKM 
observables. In this section we provide further support for this conjecture by giving some 
examples and by showing that it is possible to classify, in a procedure that is suitable for 
practical studies, the set of all mass matrices that are compatible with the given observables. 

Generally speaking, due to the use of a (as we conjecture) optimal parameterization the 
task of finding all mass matrices in NNI form is reduced to the simple problem of solving a 
quadratic equation, see (^). To begin with we note that the left-hand side of eq. (|33|) can 
be written as a scalar product 



( (m^ — m 2 c )ab , —(m 2 — m^)b, (ml — m^)a ) (h 




a 



= {m 2 t -m 2 c )ab (l, -fii/a, fjL 2 /b) (h^ la \ =0, (35) 
where the constants [i\ and \i 2 denote the ratios 

2 2 2 2 

mi — mi mi — mi 

fH = — 2 f, H2 = — 2 | = ^i-l- (36) 

mf — m A c mf — m£ 

When the condition (|35"D is written in this form and using the fact that is hermitean, we 
see at once that if (a, b) is a solution, so is (— //i/a*, ^/b*). The simultaneous substitution 

a\ f a' = -fii/a* \ , , 

b J ^ \ V = +to/V ) 

maps the circle with radius R a = ^fjxi in the complex a-plane onto itself (by relating an- 
tipodes), while in the complex 6-plane every point of the circle y^e 1 * 6 is a fixed point. 
At the same time, this substitution means interchanging the first and second columns of 
U, eq. ([[!]). Therefore, if we restrict, e.g., a = a inn0I . to the interior of the first circle, and 
calculate o(ai nne r) from ( |3"3"D as well as the mass matrices (H^'^) inner , then the solution 
pertaining to a oute r = a' — — l^i/ci* nneT and the corresponding value of fe'( a outer) yields the 
mass matrices 

(H (uUd) ) outer = Ul(H^) inneI U with 



Un 



In particular, for unprimed and primed parameters of the triangular matrices this is equiv- 
alent to 

(dinner = (Pouter Whinner = (« K 3) outcr , 
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( 7 + K 2 + K 3 ) inncr = ( 7 + K 2 + K 3 ) outcr , 

and likewise for the primed quantities: 

(a' 2 ). = (3' 2 ) , (P'n'e^). = (a'^e 1 ^) 

\ /inner V." /outer \" I /inner \ a /c 



(7' 2 + 4 2 + dinner = (V + + ^U, 



Clearly, this symmetry simplifies greatly any practical analysis. 

Before turning to the general case we illustrate our formulae by a few special cases: 
(i) If the parameter a vanishes, a = 0, eq. fl33|) gives b(a = 0) = — fai/k 2 3- For the ^-sector 
we then obtain 

- 2 2 P>2 

a = m c , p 



m l\ 


h 2 3 


2 + m t\ 


fall 


2 




^23 


2 + \ 


fax 


2 



/t 3 = , (3k 2 



\fal 


ll^23|(mf - 


™l) 


\l 


^23 


2 + 


fall 


2 



.2 = m l m \ 
7 p 2 

Similarly, for the (i-sector we find in this case 



fall 


^23 


2 + h 


33 


/? 


21 


2 _ 


- 2Re {h 12 h l3 h 2 3) 




^23 


2 


+ 


fail 


2 



v 2 





M 


1^3 1 




1^23! 


2 + 


^21 1 


2 



33 



21'*13 



/i93/i 



23 "-31 



11 



h 



23 



21 



2 = m 2 mgm 2 
a' 2 /?' 2 

This example illustrates the power and the simplicity of the reconstruction procedure: Given 
the experimental data (|14D, with Vckm given in an arbitrary, but fixed parameterization, 
the above formulae yield all entries of the triangular matrices (0), hence the mass matrices 
of the u- and <i-sectors in an NNI-basis. 

(ii) If we set 6 = 0, hence a(b = 0) = —fai/fa 2 the mass matrix in the u-sector is given by 





Ifal 


2 + m l\ 


fall 


2 




fal\ 


2 + 


fal 


2 
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For the (i-sector we obtain the following expressions 

i/2 _ u h'2 



&' 2 = h 33 , P' 





^32 


2 + h 22 \ 


hzi\ 


2 -h 


32^-13^-21 — ^-31^23^-12 






l 2 + 


^31 


2 



hi* I i#> ^32 ^31 /, , 

\hm\ l + /lai M V 



31 1 



^-31^-12 ^32^-21 

'ii i i r I ) 



32 "-31 
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As in the previous example this shows that it is possible to reconstruct the triangular 
matrices QltJ ) from the data and, from there, the squared mass matrices = T^T^ . 

The two preceding examples are degenerate cases because, by setting a (or b) equal to 
zero, hence fixing 6(0) (or a(0), respectively), the remaining two-parameter freedom is partly 
"frozen". The only freedom left over is contained in the phases ip 2 or ^3, respectively, which 
come from the phase matrix (|20D . Also, the substitution (|37D shows that two more special 
cases can be obtained where one of the parameters is sent to infinity. We also remark in 
passing that, although unrealistic in the light of the data, one can easily study the even 
more degenerate case of a and b both going to zero, for instance via 

a ^ , 6= ''i~iis ^o. 

[mi - m z u )h 2l 

In the procedure proposed by Branco et al., this limit corresponds to the case k = in 
(PI). While their analysis needs more care in this case, ours can be extrapolated smoothly 
to (a = 0, b = 0) in the way described above. Thus, there is no obstruction against choosing 
a (or b) anywhere in the complex plane. 

We now turn to the general case but keep in mind the actual values of the observables 
(quark masses and CKM angles). We first notice that with 

(m u = 5.1 MeV, m c = 1350 MeV, m t = 330000 MeV) , (38) 

(m d = 8.9 MeV, m s = 175 MeV, m b = 5600 MeV) (39) 
the first ratio fl36|) is approximately 1 while the second is very small, 
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- 1.67 x 1(T 5 , Hi = nt - 1 » 1.67 x 10~ 5 . 

It seems appropriate to expand our formulae in terms of fi 2 — A*i — 1- So, for a given value 
of a, the two solutions 61/2(0) of the quadratic equation (|33| ) are given by 

61(a) = 7 - ^ + f 3 , 2 5T a ^ + O(^) , (40) 
h 21 + a(h 22 - h n ) - a 2 h 12 

«■> = ^737^ + °w ■ < 41 > 

Whether or not this is a good approximation, in principle, depends on the range of a and on 
the matrix elements h^, hence on the experimental input. In order to estimate its quality 
it is useful to compare the product and the sum of the two approximate solutions to the 



product and the sum of the exact solutions of the quadratic equation (p3|) . Thus, denoting 
the above approximations by 6* (a), the exact solutions by 6f cact , we define 



bi(a)b 2 (a) 
6f act (a)6f act (a) 



_ = 61(a) +6 2 (a) _ ^ 

If the data are such that 5p and 5s are small, and using the fact that the ratio 61/62 is 
proportional to fi 2 , estimates for the approximate solutions are seen to be the following 

61 (°) ~1_i_Jf a ^(a) ^ i r 



exact /„\ •* ' kexact/ 



6f cact (a 

In practice, i.e. for realistic values of the experimental input, the quantities (f42"D, as well 
as the modulus of the ratio &i/6 2 are very small. Indeed, with the masses fl3"8p, flBTf), and 
with the following data for the moduli of the CKM matrix elements |6j , (assuming a positive 
value of the CP-invariant J) 

\V ud \ = 0.9752 , \V US \ = 0.2213 , |V cd | = 0.2211 , \V CS \ = 0.9744 , (43) 

the elements of the hermitean matrix = (hij) = {m 2 d + m 2 + m 2 )(kij) are found to be 

fcn = 6.144 x 10" 5 , fc 2 2 = 2.584 x 10" 3 , k 33 = 0.9973 



k 12 = (1.089 + i 2.080) x 10" 4 , k 13 = 3.352 x 10~ 3 - i 8.4890 x 10" 

k 23 = 4.062 x 10" 2 + i 6.988 x 10~ 7 . 
The quantity 5p is easily seen to be 

dp [a) = : 7—^2- 

ak 13 - k 23 
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As the numerical values of the elements of are such that I/C23I/I&13I ~ 12 this function 
is regular over the interior of the circle with radius R a = ^f]i{ and may thus be estimated 
by means of standard techniques of function theory. We find 

\5 P \ < 1.82 x 1(T 5 . (44) 

Due to the symmetry ( |37|) this domain is sufficient to cover all NNI solutions. Estimating 
Ss(a) is a bit more complicated because, as a function of a, it has two poles in that same 
domain, very close to each other. Excluding a small circle around these poles one finds 
typically 

\6 S \< 0.023. (45) 
Also the ratio &1/&2 as obtained from (|40|) and ( fiTf ) is estimated as follows 



< 0.032 



Given the experimental values (^), d39|) , and (|43|) , 5s, eq. ([4~5"D is the dominant uncertainty. 



Thus, in this framework the expressions fl40|) and (fflj) are excellent approximations in the 
interior of the circle with radius R a except in a small neighbourhood of the two poles of 
61(a), eq. (fK]). Note that the approximations are continuous in the parameter a which is to 
say that the reconstructed mass matrices and depend on the remaining freedom 
in a continuous manner. This is particularly relevant when studying the dependence of the 
mass matrices on the parameter a and when comparing to textures obtained from specific 
physical assumptions. 

We illustrate the method by means of two examples in Figs. 1 and 2. These figures show 
the parameters a 2 and (3 2 as functions of the complex parameter a and for the two approxi- 
mate solutions 61/2(0), eqs. (|4*0"| ) and ([H]), with a chosen from the interior of the circle with 
radius R a . The figures show clearly the smallness of the neighbourhood of the two poles 
where the approximation breaks down. Note that if in that region one wishes to use the 
exact solutions fef^ * care must be taken in insuring continuity when the signs of square 
roots are chosen. 



V. SOME EXAMPLES AND CONCLUSIONS 

In a detailed numerical study M we have verified that the procedure that we are proposing, 
from a practical point of view, is manageable and transparent, and that all dependencies of 
the NNI parameters of the triangular matrices can be illustrated in a simple manner. 

Assumptions about specific textures of the mass matrices obtained on the basis of some 
physical conjecture, may or may not be compatible with the data. The parameterization 
that we are proposing in this work is particularly well suited for testing the consistency of 
any such model in a simple and transparent manner. We illustrate this statement by two 
examples taken from the literature. Suppose M in an NNI basis the w-sector is constrained 
to be, in addition to the NNI condition, 
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/ o n o \ 

M mhei = I r i r 2 , with rx,r 2 ,r 3 real. 
\ r 2 r 3 / 

The corresponding hermitean, squared form is 

(r\ rir 2 \ 
r? + rf r 2 r 3 
^1^2 ^2^3 r l + r l J 

with 7-3 = + m 2 c + — 2{r\ + r 2 ). Comparing this to = T^T^ we see that 

-2 2 52 2,2"- a* 

a = r 1 , fj = r 1 +r 2 , ozk 3 = r x r 2 , pn 2 = r 2 r 3 . 

Thus, by eqs. (|27|) , (p8|) the moduli of a and b are fixed and the constraint ( |33D is reduced to 
an equation for the pase factors e"^ a and e 1 ^ 6 . It is then easy to decide whether or not this 
equation has a solution and, thereby, whether or not the ansatz of the model is compatible 
with the data. 

The second example || makes again use of an NNI basis but now constrains the <i-sector 
further by assuming 

/ Sl \ 
M modei= U «2 . witn Sj, s 2 , s 3 real. 
V s 3 s 3 J 

Like in the previous example the remaining freedom is reduced to two phase factors which 
must obey the constraint (|33p. As the latter contains the input data, i.e. quark masses and 
CKM mixing angles, it is not clear a priori that the model is admissible. We note in passing 
that according to |§ both models, within the experimental error bars, can indeed be used 
to parameterize the data. This is checked in our framework by confirming that eq. ( j33|) has 
solutions of modulus 1. 

The point we wish to make by quoting these examples is the following: while in general 
it is difficult to test the compatibility of a specific model ansatz with the data (wihin their 
experimental error bars), the model may always be transformed to an NNI-basis. By con- 
verting it to our general form in terms of the parameters a and b, its test in the light of the 
data is reduced to checking the simple quadratic equation 

In summary, we found a new parameterization of squared mass matrices in terms of the 
experimental input (eigenvalues and mixing observables) and one complex parameter that 
allows to sample the space of solutions in an analytical and transparent manner. Indeed, 
from the input: quark masses, matrix elements hij as obtained from the CKM data, eq. ([T1D, 
and a choice of a (from which b(a) is obtained via (|33|), or vice versa), the equations given 
in appendix B directly yield the mass matrices fllpp . Thus, by varying the parameter a over 
the circle with radius R a in the complex plane, and using the symmetry fl37D we scan the 
space of all admissible mass matrices, up to unobservable changes of basis. 

We conjecture that this procedure of reconstructing all mass matrices which are compatible 
with the data up to (unobservable) changes of bases, is optimal. We obtained this result 
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by combining the idea of using general NNI-bases |T| with the polar decomposition theorem 
that allows to restrict the general analysis to triangular matrices @], ||. The formulae 
that we obtained are sufficiently simple to handle so that they may be implemented in a 
reconstruction routine that also takes account of the experimental error bars. Alternatively, 
as demonstrated by the examples we gave, our method allows for a quick test of compatibility 
with the data for any assumed texture in the mass matrices. 

Finally, with our knowledge of neutrino oscillations and of the corresponding mixing 
matrix increasing, it will eventually be possible to perform the analogous analysis of the 
leptonic mass matrices in the standard model. 



APPENDIX A: 



This appendix gives some intermediate results which are skipped in the main text of 
section III. We begin with the expressions for k 2 and k% in terms of a and (3. 
Inserting k\ according to (^3) into (|22| ) and making use of leads after a straightforward 
calculation to 



h ■> 



In a similar way we also get: 



k\ 



(a 2 — m 2 u ){m 2 — 6?)(m\ — a 2 ) 
a 2 2 - a 2 ) 



2 -m 2 u )0 2 - m 2 c )(m 2 -p 2 ) 



(3 2 ((3 2 - a 2 ) 

Defining U =: PV, where U is given by and V = (v^) we thus obtain: 
/ 



V i3 



2 - 


ml){m 2 c - 


a 2 )(m 2 - 


-d 2 ) \ 






( _ 


(m 2 - 


-ml){m%- 


-m 2 )(f3 2 


-a 2 ) 








(a 2 - 


ml)W 2 - 


m 2 )(m 2 - 


-a 2 ) 




Vi2 = 


± 


K- 


-m 2 )(m 2 - 


-m 2 )0 2 


-a 2 ) 


) 


(a 2 - 


m 2 )(m 2 - 


-a 2 )(m 2 - 












-m 2 ){m 2 - 


-m 2 )0 2 


-a 2 ) J 






V 


in 


-m 2 )0 2 - 


-ml) 




(m 2 


-a 2 )0 2 - 


ml) 




-m 2 )(m 2 


-m 2 ) ' 






-m 2 )(m 2 - 


-ml) 



(a 2 - 


-m 2 u )0 2 - 


-m 2 )(m 2 - 


-/3 2 ) \ 


(m 2 - 


-m 2 )(m 2 


-m 2 )(/3 2 


-a 2 ) 


2 - 


-m 2 )(m 2 - 


-& 2 )(m 2 - 


-P 2 ) 


(mf- 


-m 2 )(m 2 


-m 2 )0 2 


-a 2 ) 


(/3 2 - 


-m 2 )0 2 - 


-m 2 )(m 2 - 


-a 2 ) 


(m 2 - 


-m 2 )(m 2 


-m 2 )0 2 


-a 2 ) J 








\ T 



(m 2 -a 2 )(m 2 -/3 2 ) 
(m 2 -m 2 )(m 2 -m 2 ) 



In (|A3|) the upper signs for v 2 ± and v 2 2 refer to the case 

m„ < a 2 < m 2 and m 2 c < p < m 2 (case I) 

whereas the lower signs pertain to the case 

ml < a 2 < m 2 and m 2 u < j3 2 <m 2 (case II) 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



In fact, this distinction of two cases easily follows from the positivity of n 2 , i = 1,2,3, see 
At the same time, this argument shows that these two cases exhaust all possibilities. 
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APPENDIX B: 



In this appendix we quote the results for the parameters of the mass matrices in NNI 
form in terms of the complex variables a and b (|26]). For the up-sector the expressions in 
question are obtained by means of ([18]) and inserting U according to (0), or, equivalently, 
by using |§) and Q - (T£3]): 



a 



Ni 



(|a| 2 |6|XK - m 2 c ) 2 + \a\ 2 m 2 t (m 2 c - m 2 u ) 2 + \b\ 2 m 2 c (m 2 - m 2 u f) 



/3 2 
(3k 2 



(ml + \a\ 2 m 2 c + \b\ 2 m 2 ) 

N2 
N 2 



&k 3 = j^\a\\b\(mt -m c )(m t -m u )(m 



2 2^ 
c ~ m u, 



(Bl) 
(B2) 

(B3) 
(B4) 



In addition, 7 2 is determined from (0), viz. 

„ 2 (m u m c m t ) 2 
7 = 

a 2 (3 2 

The results for the down-sector follow from 

#W = UH^U ] , 
where is given in fll4]) . Setting = (hij) we thus get: 

'Si 



(B5) 



a 



j^{(mi - m 2 c Y\ab\ 2 h n + (m 2 t - m 2 u ) 2 \b\ 2 h 22 + (m 2 c - m 2 u ) 2 \a\ 2 h 33 

2{m 2 t - m 2 u ){m 2 - m 2 )|6| 2 Re {ah 12 ) + 2(m 2 - m 2 u ){m 2 - m 2 )|a| 2 Re {bh 13 ) 
2(m 2 c -m 2 u )(m 2 -m 2 u )Re(abh 2 , 



(B6) 



/3' 2 



h n + \a\ 2 h 22 + 



\b\ 2 h 33 + 2 (Re (ah 12 ) + Re (bh 13 ) + Re (abh 23 )) } (B7) 



1 



y-2 y., t O lVc - m lf + l & ! 2 K - 2 ](^H + fl/i 12 + ^13) 

[|6| 2 (m 2 - m 2 ) - (m 2 - m 2 )]a*(/i 21 + bh 22 + 6/i 23 ) 
+ [|a| 2 (m 2 - m 2 ) + (m 2 - m 2 u )}b*(h 31 + a/i 32 + bh 33 ) 
a(h 3 i + a/i 32 + bh 33 ) - b(h 2 i + a/i 22 + bh 23 ) N 2 



ab(m 2 — m 



N 2 



(B8) 
(B9) 



1 



a K 3 e 



N x N 2 \ab\ 



[|a| 2 K 2 - ml) + \b\ 2 (rn 2 - m 2 )] [-\ab\ 2 {m 2 - m 2 )/^ 



+ a|&| 2 (m 2 - m 2 u )h l2 - %| 2 (m 2 - m 2 u )h 13 ] 
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+ [\b\ 2 (m 2 t - ml) - {m 2 c - m 2 u )][-a*\ab\ 2 (m 2 - m 2 c )h 2l 

+ \ab\ 2 (m 2 - m 2 u )h 22 - a*b\a\ 2 (m 2 c - m 2 u )h 22 \ 

+ [\a\ 2 (m 2 - m\) + (m 2 - m 2 )] [b*\ab\ 2 (m 2 - m 2 c )h 31 

- ab*\b\ 2 (m 2 - m 2 u )h 32 + \ab\ 2 {m 2 c - m 2 )/i 33 ]} ■ (BIO) 

f 2 

The second, alternative form ( p9|) of 



_ (m d m s m b ) 
a' 2 (3' 2 

(|B8|) is obtained by a few lines of calculations. 
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FIGURES 





FIG. 2. Part a: the parameter p 1 as a function of a (real and imaginary part) for the first 
solution (|40|); Part b: same parameter for the second solution (fll|) 
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